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The knowledge of the nuclear level density is necessary for understanding various reactions in¬ 
cluding those in the stellar environment. Usually the combinatorics of Fermi-gas plus pairing is 
used for finding the level density. Recently a practical algorithm avoiding diagonalization of huge 
matrices was developed for calculating the density of many-body nuclear energy levels with certain 
quantum numbers for a full shell-model Hamiltonian. The underlying physics is that of quantum 
chaos and intrinsic thermalization in a closed system of interacting particles. We briefly explain this 
algorithm and, when possible, demonstrate the agreement of the results with those derived from 
exact diagonalization. The resulting level density is much smoother than that coming from the 
conventional mean-field combinatorics. We study the role of various components of residual inter¬ 
actions in the process of thermalization, stressing the influence of incoherent collision-like processes. 
The shell-model results for the traditionally used parameters are also compared with standard phe¬ 
nomenological approaches. 

PACS numbers: 21.10.Ma, 21.10.Hw, 21.60.Cs 


I. INTRODUCTION 

The knowledge of the level density is an important ele¬ 
ment in understanding the behavior of a quantum many- 
body system of interacting particles in various physical 
processes. In nuclear physics, this knowledge is necessary 
for the description of numerous reactions, including those 
of astrophysical or technological interest. The cross sec¬ 
tions can be very sensitive to the level density that typi¬ 
cally grows exponentially as a function of the excitation 
energy and the number of constituents. In turn, the the¬ 
oretically predicted level density is sensitive to the statis¬ 
tics of particles, their specific interactions and available 
orbital space that, in realistic computation, usually has 
to be truncated. Apart from that, the level density in a 
finite self-bound system, such as the atomic nucleus, can 
be different for the available classes of eigenstates charac¬ 
terized by different quantum numbers of exact constants 
of motion (in nuclei total spin J , parity II if we neglect 
weak interaction, and isospin T if we neglect interactions 
violating charge-independence). 

Below we consider a problem of practical microscopic 
calculation of the level density for a nucleus described 
by a Hamiltonian of the shell-model type. In this frame¬ 
work it does not matter if the Hamiltonian is derived 
from a more fundamental approach or fit phenomenolog¬ 
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ically with the use of experimental data. We assume that 
the Hamiltonian describes the low-lying energy spectra, 
transition rates and other observables known from the 
experiment reasonably well. Then the task is to predict 
the level density of the system at higher excitation en¬ 
ergy, in the region beyond direct measurements resolving 
individual quantum states. Practically, the relevant fac¬ 
tual milestones, apart from the low-lying spectroscopy, 
are the regions of isolated neutron resonances near neu¬ 
tron separation energy and the results of the Oslo method 
and related experimental approaches M. Of course, 
any practical shell-model Hamiltonian loses its validity 
outside of the truncated orbital space where this Hamil¬ 
tonian was expected to work. At some excitation en¬ 
ergy, the states of the system come from the particle or¬ 
bitals not included in the model. However, with available 
computational means, the space of validity of the model 
Hamiltonian can be broad enough, in particular including 
the excited states involved, for example, in astrophysical 
reactions at a typical stellar temperature. We can also 
hope to use the microscopic results for the nuclei far from 
stability, where the level density is usually predicted by 
pure phenomenology [B|. 

Another physical limitation arises from the obvious 
fact that the states involved in the reactions belong to 
the continuum, while the standard shell-model calcula¬ 
tions work in the discrete spectrum. Instead of discrete 
states here one has to deal with resonances seen in var¬ 
ious reactions. Then the whole definition of the level 
density becomes questionable and, strictly speaking, one 
has to move to the complex plane of resonances. How- 
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ever, the traditional approach is still meaningful if the 
typical widths of the involved states are small compared 
to the spacings between the states with the same quan¬ 
tum numbers. In what follows we limit ourselves by this 
situation. 

Neglecting the continuum effects, the trivial solution 
for the level density generated by a certain Hamiltonian 
would be a full diagonalization of the Hamiltonian matrix 
in an appropriate orbital basis. However, this is practi¬ 
cally only possible in sharply truncated orbital spaces, 
which might be only sufficient for relatively light nuclei, 
like in the sd-shcll Q. In many realistic cases of cur¬ 
rent interest, the dimensions of corresponding matrices, 
even in subspaces with given quantum numbers, are pro¬ 
hibitively large. Moreover, such a diagonalization is any¬ 
way superfluous because we do not need full information 
on every excited state in spectral regions of high level 
density. The level density is, by construction, a statisti¬ 
cal notion. 

In this situation, we are looking for the statistical so¬ 
lution of the problem. The physical justification of such 
an approach lies in the fact that, at small level spac¬ 
ings, the stationary nuclear states are extremely compli¬ 
cated superpositions of simple determinantal states with 
integer occupation numbers of definite orbitals. Gradu¬ 
ally switching on inter-particle interactions and going in 
this process through multiple avoided crossings of vari¬ 
ous configurations we come to chaotic states [6, 7j with 
observable properties smoothly changing along the spec¬ 
trum. Therefore our problem reduces to finding a realis¬ 
tic way to describe this smooth evolution. 

This purpose can be reached using the methods of 
statistical nuclear spectroscopy [8l l 1 (j| . Already in the 
framework of a single partition (a certain configuration 
of independent particles occupying the mean field levels), 
the level density after including the particle interaction 
rapidly goes to the Gaussian limit with the increasing 
particle number [llj . This is some kind of manifestation 
of the central limit theorem. The many-level, and there¬ 
fore many-partition, generalization should give a reliable 
image of the total level density for an accepted orbital 
scheme. This has to be done for each class of global 
quantum numbers. This direction of theoretical search 
has a long history. We would like especially mention the 
works in the direction of statistical spectroscopy applied 
to shell-model Hamiltonians, see for example [H, fl3j . 
After several preliminary publications, our successful al¬ 
gorithm was constructed [14| and opened for public use 
fl5| . The results of implementing this algorithm for the 
level density in sectors with given values of global con¬ 
stants of motion are practically identical to those from 
the full diagonalization when the latter is possible. For 
well-tested shell-model Hamiltonians, the results are in 
good agreement with available experimental data. 

For many years, starting with the classical work by 
Bethe [H], the nuclear level density was estimated us¬ 
ing the combinatorics based on the ideas of a Fermi gas. 
The influential review of earlier approaches of this type 


was give n by Ericson 0 ,thc later derivations can be 
found in (l8lT2lj| . see also [221. The recent achievements 
in this direction [23| - |25[ include the pairing correlations 
considered as a part of the self-consistent mean field 
in the framework of the BCS theory or Hartree-Fock- 
Bogoliubov variational ansatz. The shell-model Monte- 
Carlo methods [26l - [28| . being very demanding computa¬ 
tionally, work relatively well at least with some parts of 
the full shell-model interaction but require the projection 
to the correct values of spin and parity. 

The chaotization of the dynamics mentioned above 
leads to the possibility of describing the physics of excited 
states at high level density in terms of statistical thermo¬ 
dynamics including temperature, entropy etc. This was 
understood in application to nuclear reactions from the 
early times of nuclear physics [Tol |29l H(j . The detailed 
analysis of atomic [3 lL l32j and nuclear [g], [33[ chaotic 
states supported an old idea Q of thermalization in a 
closed system driven by the interactions between the con¬ 
stituents, with no heat bath: the average over a generic 
chaotic wave function in a chaotic region is equivalent 
to the average over a standard equilibrium thermal en¬ 
semble [35f. Currently this idea, sometimes called the 
“eigenfunction thermalization hypothesis”, is extensively 
discussed in the many-body physics community [36| . One 
of the purposes of the current publication is in compar¬ 
ison of the exact shell-model nuclear level density with 
phenomenological ideas based on the Fermi-gas picture at 
certain temperature. We look at these ideas and based 
on them equations from the viewpoint of our numeri¬ 
cal results. Our attention will be mostly concentrated 
on the usually cited empirical parameters of the level 
density and their energy and spin dependence. Another 
point of interest is in the role of various components of 
the shell-model interactions in the formation of the level 
density. One important result is that the consideration 
of the mean field, even with addition of the BCS-type 
pairing, is not sufficient. Incoherent components of the 
interaction in a finite many-body system, as a rule ne¬ 
glected in the mean-field combinatorics, play a significant 
role smoothing the energy behavior of the level density. 

In what follows we briefly explain the method and give 
the examples of practical calculations. The results will 
be compared with what would follow from traditional 
phenomenological models. 

II. MOMENTS METHOD 

We consider a finite system of interacting fermions de¬ 
scribed by the standard Hamiltonian 

H = ^2 + — ^2 Tl2;340i(4 a 3 a 4 (1) 

1 ^ 1234 

that contains the mean-field part with effective single¬ 
particle energies e\ and the antisymmetrized two-body 
interaction. The generalized numerical subscripts com¬ 
bine all quantum numbers of single-particle orbitals. In 
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this form the method can be applied to nuclei, atomic 
or molecular electrons, atoms in traps etc., any system 
where the residual interaction is sufficiently strong to 
produce complicated eigenstates. The three-body forces 
can be included in the same way although the computa¬ 
tions become more cumbersome. Using a phenomenolog¬ 
ical shell-model Hamiltonian we assume that many-body 
forces, at least partly, are included in fitted matrix ele¬ 
ments. The quality of the Hamiltonian is checked by the 
explicit applications to individual low-lying states and 
comparison with available experimental information. 

In a finite self-bound system, such as the atomic nu¬ 
cleus, the total angular momentum (nuclear spin) is ex¬ 
actly conserved supplying good global quantum numbers 
J and M. Therefore it is convenient from the very be¬ 
ginning to use a spherically symmetric basis of single¬ 
particle orbitals | jm) which define, along with the orbital 
momentum £, main quantum number is, and isospin r, 
the quantum numbers combined in Eq. (ID into a unified 
numerical subscript. If the orbital space is sufficiently 
broad, this spherical shell model can describe intrinsic 
deformation without violating rotational symmetry [37f. 
In practice, the operators in the Hamiltonian can be com¬ 
bined into pairs, such as (a^apLA with the angular mo¬ 
mentum quantum numbers of the pair, LA, fixed through 
the Clebsch-Gordan coefficients. In the same way, for the 
isospin-invariant forces, the isospin of the pair can be also 
fixed, and then the interaction V in a restricted orbital 
space is defined through a finite number of pairwise ma¬ 
trix elements. If more convenient for computations, one 
can as well use the M-scheme without vector coupling. 

The practical algorithm of calculations in the M- 
scheme follows from individual configurations, p, ( par¬ 
titions ) which are possible distributions (ni, ri 2 ,...) of 
available particles, X^i n i = N, over single-particle or¬ 
bitals (here the index 1 does not include the projection 
m). The many-body states |a) possible for each parti¬ 
tion form a subspace where a combines total quantum 
numbers N, M, T 3 and parity. It is convenient 14[ to use 
the proton-neutron formalism. 

Let D ap be the dimension of the class of states with 
global quantum numbers a built on the partition p. As 
shown in statistical spectroscopy SE3 and confirmed in 
many examples by the exact shell-model diagonalization, 
the density of states for a given partition is close to the 
Gaussian. Of course, this is the main assumption based 
on a rich experience with the features of quantum chaos 
in mesoscopic systems. The characteristics of the Gaus¬ 
sian are defined by the moments (traces) of the actual 
Hamiltonian. The centroid is just the mean energy value 
for a given partition, 

E ap = {H) ap = Tr ( “ p) M. (2) 

The dispersion of the Gaussian, o ap , is the second mo¬ 
ment, 

a 2 ap = (H% p - E 2 ap = -L Tr ^ H 2 - E 2 ap . (3) 

-E'ap 


It is important to stress that the calculation of these 
traces does not require the diagonalization of large-scale 
matrices. The first moment m is the diagonal matrix 
element of the Hamiltonian averaged over the partition, 
while the second moment © is the sum of squared off- 
diagonal elements along one line of the Hamiltonian ma¬ 
trix, again averaged over the lines corresponding to the 
partition. It is known that the dispersion for each ba¬ 
sis state very weakly fluctuates within a partition jf| HU 
even prior to the next averaging. The second moment 
includes all interactions coupling the partitions. If the 
traces were calculated in the M-scheme, we obtain the 
density of states counting all M- -degenerate states within 
the nrultiplets. To obtain the level density for given spin 
J, we have in a standard way to find the difference of 
traces for M = J and M = J + 1. 

The total level density is given by summing the con¬ 
tributions of partitions using the constructed Gaussians 
G ap {E) with their centroids © and widths ©, 

p{E\a) =^D ap G ap (£). (4) 

V 

As was understood earlier 3S!], it is better to use finite 
range Gaussians, 


G ap {E ) — G(E — E ap + -Eg.s.; crap), 


where 


G(x;ct) =C 


g-* 2 / 2 ^ \ x \< w 

0, |x| > TjCT 


( 5 ) 

( 6 ) 


Here E s , s . is the ground state energy to be defined sep¬ 
arately, C is the normalizing factor, f dxG(x,o) = 1, 
and the cut ting finite-range parameter 77 has to be found 
empirically ©j; its actual value rj « 2.8 agrees with the 
analysis of the shape of typical nuclear strength functions 


Following the recipe of Ref. [42| we make an important 
addition to the algorithm. In many realistic applications, 
the orbital space contains several shells labelled by the 
harmonic oscillator quanta N . In this case, the standard 
formulation of the shell model includes cross-shell tran¬ 
sitions with unphysical excitations of the center-of-mass. 
These spurious states are to be excluded from the level 
density. In some versions of the shell model, these states 
are artificially shifted to high energies. Here, the sub¬ 
traction of ghost states is accomplished by renormalizing 
the contaminated level density p(E,J;AT) through the 
recurrence relations. For example, while the Af = 0 case, 
which will be called the p°(E, J; 0) approximation, is free 
of admixtures, the pure level density p° at the next step 
(no admixtures of the single center-of-mass excitation) is 
found as 


J +1 

P°(E, J; 1) = p{E, J; 1) — ]T p(E,f ;0). (7) 

1 | 
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FIG. 1: Nuclear level densities for 28 Si, II = +1, various spins; shell model for the sd-shell and USDB two-body interaction 
(solid curves) vs. moments method (dashed curves); finite-range parameter rj = 2.8. 


Here the sum goes over the intermediate angular mo¬ 
menta J' from | J — 11 to J + 1 since the center-of-mass 
operator is equivalent to a vector. If higher admixtures 
M > 1 are present, the recurrence relation has to include 
a corresponding number of steps back. This makes the 
calculation of the trace of H 2 with various intermediate 
states slightly more involved [40], |43j . 


III. EXAMPLES OF LEVEL DENSITY 

A. Comparison with the exact solution of the shell 
model 

The first natural check of the approach is in compar¬ 
ison of the resulting level density with the picture aris¬ 
ing from the full shell-model diagonalization in the cases 
where such a diagonalization is technically plausible. The 
sd-shell model for a long time is known as the best ex¬ 
ample of exact diagonalization. The model is completely 
fixed by the effective single-particle energies d$/ 2 ,Si/ 2 , 
and ^ 3 / 2 , and 63 phenomenologically fitted matrix el¬ 
ements of two-body interaction. The model works ex¬ 
tremely well for practically all observables of sd-nuclei 
and not only for the lowest states. For ex amp le, both, 
the experiment (44} , and the sd-shell model [42| . indicate 


the existence of ten stationary states with J u = 0 + up 
to excitation energy of 15 MeV in 28 Si, therefore provid¬ 
ing the same average level density at least at not very 
high energy; the mean level spacing between those 0 + 
levels is 0.95 MeV in the experiment and 1.02 MeV in 
the shell-model calculation. 

Fig. [T] illustrates the results of the moments method 
for calculating the level density in 28 Si, a typical object 
of the sd-shell model applications that has served long 
ago as a testing ground for quantum chaos Q. We see 
that the level density for different classes of states, here 
0+,l + ,2+,and 3+, is always a smooth curve of the Gaus¬ 
sian type. Of course, as the calculations have been done 
in the restricted orbital space, the real physical result 
that can be juxtaposed to the experimental data and used 
for the reaction calculations always corresponds only to 
the left hand side of the full graph and to the excitation 
energy below the centroid maximum. 

Various specific quantum numbers (total spin and par¬ 
ity) produce the level density of the same qualitative be¬ 
havior, with the integral corresponding to exact multi¬ 
plicities of states with a given set of quantum numbers 
in a fixed orbital space. All examples look the same. The 
agreement with the exact shell-model diagonalization is 
almost perfect, with slightly more visible fluctuations for 
the class J n = 0 + that has a smaller total dimension. In 








FIG. 2: (Color online.) Comparison of the nuclear level densities of 20 Ne in lftcj-approximation, II = —1. The densities are 
calculated within the shell model approach (stair-case curves) and with the moments method (solid curves). The black curves 
on the left and right graphs (both are higher at low energies) correspond to the densities with spurious states included; the 
blue curves (that are lower at low energies) correspond to the densities without spurious states; the red (dotted) curves present 
the spurious nuclear level densities. 




Excitation energy (MeV) Excitation energy (MeV) 

FIG. 3: (Color online.) Comparison of nuclear level densities calculated with the shell model for 22 Mg in lhcj-approximation, 
II = —1, and the moments method. The coloring scheme is the same as in Fig. [2] The small peaks around 140 MeV show the 
spurious densities. The inserts present the low-energy part. 


all cases we see a small deviation near the centroid which 
supposedly can be eliminated by taking into account the 
fourth moment of the Hamiltonian (but it makes no sense 
to go for such complicated and time-consuming calcula¬ 
tions to improve the results in the region outside the 
physically relevant area). 

The smoothness and Gaussian behavior of results in all 
cases confirm the possible thermodynamic interpretation 
in terms of entropy S(E) (mean logarithm of the level 
density) and temperature, dS/dE = 1/T(E). Formally, 
the centroid of the level density for the finite orbital space 
corresponds to infinite temperature and the right half 
of the curve to negative temperatures. The full shell- 
model analysis of the wave functions Q has found that 
the same effective temperature can be extracted by the 
single-particle thermometer using the occupation num¬ 


bers of available spherical orbitals for individual station¬ 
ary many-body states. The interaction of the quasipar¬ 
ticles in the self-consistent mean field acts as the heat 
bath, and the chaotic mixing of the eigenstates leads to 
thermalization even in such a small Fermi system. 


B. Elimination of spurious states 

As mentioned above, the full shell-model diagonaliza- 
tion in the cases with the presence of transitions between 
the orbitals of opposite parity (excitations across the os¬ 
cillator shells) brings in ghost states related to the center- 
of-mass motion rather than to intrinsic excitations. We 
explained above the recurrent techniques used for elimi¬ 
nating these spurious states and obtaining the pure level 





























FIG. 4: Comparison of the experimental nuclear level densities (stair lines) and the densities calculated with the moments 
method (straight lines) for 28 Si, all J, positive (the left graph) and negative (the right graph) parity. 




FIG. 5: Comparison of the experimental nuclear level densities (stair lines) and the densities calculated with the moments 
method (straight lines) for 26 Al, all J, positive (the left graph) and negative (the right graph) parity. 


density. The recipe frequently used in the shell model is 
the brute-force shift of the undesired states to high en¬ 
ergy by adding to the Hamiltonian under diagonalization 
a Lawson term J4f| that in the harmonic oscillator field 
of frequency oj looks as (/3 > 0) 


H' = 0 


A 

Tloj 


( 8 ) 


As was shown long ago Q, this recipe indeed generates 
a new branch of eigenstates shifted to high energy (by 
about ~ PMA) but having essentially the same com¬ 
plexity (measured by the information entropy) as their 
predecessors without spurious admixtures. 

The separation of unphysical spurious states according 
to the recurrence relations 0 works well. The standard 
shell-model shift of these states to higher energy pro¬ 
duces exactly the same group of states as our procedure. 
Fig. [2] shows the level density calculated for 20 Ne in the 
broad space of s + p + sd + pf orbitals in the Vhuj ap¬ 


proximation (only one-step excitation of negative parity 
states that include non-physical admixtures). The full 
(“blind”) calculation compared to that treated with re¬ 
currence relations 0 contains the excess shown by the 
small (red dotted) Gaussian-type curves at relatively low 
energy. We see that this curve of difference exactly coin¬ 
cides with the result of the shell-model shift according to 
Eq. 0. This means that the method of recurrence rela¬ 
tions correctly removes the unphysical excitations. The 
same conclusion can be made for 22 Mg, Fig. 0 where the 
insets zoom out the region of low energy. 


C. Comparison to experimental data 

As stated earlier, the density of low levels with J n = 
0+ in the favorite nucleus 28 Si of the sd shell model is in 
good agreement with the data. Figs. 0] and [5] compare the 
level density calculated with the moments method (solid 
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FIG. 6: Nuclear level densities for 52 Fe and 52 Cr, II = +1, different spins: shell model calculation (solid curves) vs. moments 
method (dashed curves) and Hartree-Fock + BCS (dotted corves). Finite-range parameter r] = 2.8, p/-shell, GXPF1A 
two-body interaction. 


lines) with the experimental results (stair-cases). The 
solid stair-cases present an “optimistic” attitude when 
all experimental levels of uncertain parity were counted. 
Opposite to that, the dashed lines exclude the levels 
whose parity is given by the experiment only tentatively. 
The left graph of Fig. [4] shows the summed level den¬ 
sity of all positive parity states in 28 Si calculated for the 
s + p + sd + pf space up to 14 MeV excitation energy. 
The same quality of comparison can be seen for the pos¬ 
itive parity states in 26 Al, the left graph of Fig. 0 the 
odd-odd nucleus with the relatively well measured energy 
spectrum. 

For the negative parity states, the right parts of Figs. 
[]] and [5J the shell model predicts more levels than until 
now have been found experimentally. This can be both 
due to the imperfection of the shell-model Hamiltonian 
and because of incompleteness of the data. 


IV. SHELL-MODEL PREDICTIONS AND 
MEAN-FIELD COMBINATORICS 

The widely used standard road to the nuclear level 
density is going through the mean-field representation of 
the nuclear dynamics. This traditional approach is based 


on the classical idea EM3 of the Fermi-gas where the 
excited levels result from combinations of many particle- 
hole excitations. Practically, the combinatorics is used of 
single-particle excitations from the fully occupied Fermi 
surface identified with the ground state population of the 
lowest individual orbitals. At low excitation energy, a 
renormalization of the level density is related to the gap 
due to the Cooper pairing. In complex nuclei, the low- 
energy levels observed inside the gap can be interpreted 
as collective excitations, vibrational or/and rotational. 
As the collective phenomena of these types correspond 
typically to the slow self-consistent motion of many par¬ 
ticles, it is natural to expect that such coherent combina¬ 
tions of single-particle excitations partly compensate the 
deficit of levels at low energy due to the pairing gaps and 
give rise to the so-called collective enhancement of the 
level density [22], |46| in comparison to the single-particle 
combinatorics of independent particles and holes. Mod¬ 
ern refined approaches of this class account in various 
forms for the pairing phenomenon that changes the exci¬ 
tation spectrum, especially in even-even nuclei |23l - l25j . 

In the spirit of the mean-field combinatorics, one has to 
expect the corresponding suppression of level density at 
higher excitation energy (damping of collective enhance¬ 
ment); the level density is just redistributed. When the 
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FIG. 7: Level densities for 64 Ge, spin J — 0 and II = +1. The solid curve presents the calculation in the pf shell with the 
GXPF1A interaction, the dashed curve corresponds to the calculation in the larger model space with the level gg/ 2 added, 
the dotted curve on the left graph presents the results obtained using the HFB single-particle energies and the combinatorial 
method. 




FIG. 8: (Color online.) Level densities for 28 Si, sd model space. Different curves correspond to different scale factors: 
k\ = k 2 = {0.1, 0.2, 0.3, 0.5,1.0} when the pairing and non-pairing parts of the interaction scale similarly. The left graph 
corresponds to the total density with all J included, while the right graph describes the evolution of the J = 0 density. 


general level density grows, the vibrational modes be¬ 
come strongly mixed with simpler excitations of the two- 
quasiparticle and more complicated structure, as known 
very well from the widths of the giant resonances. With 
smoothing shell gaps, it is harder to distinguish between 
rotational and intrinsic motion. Recent experiments in 
nuclei, where low-lying collective excitations are well 
known, did not find phenomena of collective enhance¬ 
ment and its fade-out H3] ■ 

In not too heavy nuclei, the quasiparticle combi¬ 
natorics (on the base of the BCS or Hartree-Fock- 
Bogoliubov pairing description) reveals step-wise effects 
of subshell occupation and pair breaking. This leads, at 
relatively low energy, to the irregular picture of the level 
density that clearly reflects these steps. The shell-model 
Hamiltonians, as a rule, contain all interaction matrix el¬ 


ements allowed by the selection rules. One of the main 
conclusions of the full shell-model calculation is that the 
presence of all interactions is significantly smoothing the 
whole picture so that it is hard to see the traces of in¬ 
dividual families which could be still recognized only by 
the special observables and selection rules for the indi¬ 
vidual transitions. We can recall that our algorithm still 
starts with the partitions formed by independent parti¬ 
cles which then overlap and lose their boundaries. 

The shell-model Hamiltonian contains all pairing ma¬ 
trix elements (and not in the simplified form with con¬ 
stant matrix elements) as well as the interaction pro¬ 
cesses responsible for multipole-multipole forces and de¬ 
formation. Therefore all collective effects mentioned 
above are fully taken into account if the orbital space 
is sufficiently broad. Fig. [G] shows the comparison of 






























FIG. 9: (Color online.) Level densities for 28 Si, sd model space. The non-pairing interaction is always off: &2 = 0.1, while the 
pairing interaction scales, fci = {0.1, 0.2, 0.3, 0.5, 1.0, 1.5}. The left graph corresponds to the total density with all J included, 
while the right graph describes the evolution of the J = 0 density. 




FIG. 10: (Color online.) Level densities for 28 Si, J = 0, 
sd model space. The black curve presents fci = 1.0, fe = 
0 .1, then the remaining parts of the interaction are in¬ 
creased together with up to the red curve (fci = = 

1 .0) that shows the density for the realistic interaction. 


FIG. 11: (Color online.) Level densities for 52 Fe, all J, pf 
model space. This figure and its color scheme are similar 
to the left panel on Fig. [8] 


the level densities for the states 0 + and 1 + in the nu¬ 
clei 52 Fe and 52 Cr of the p/-shell. The thin dotted lines 
give the level density found with the mean-field combi¬ 
natorics built on the Hartree-Fock mean field and BCS 
pairing description. All irregularities of the level den¬ 
sity found through the mean-field combinatorics are com¬ 
pletely smoothed in the full moments calculation. This 
is a typical result encountered in all examples. Again 
we see that the method under discussion produces the 
level density practically identical to the full shell-model 
diagonalization when the latter is possible. 

Fig. [7] illustrates the influence of the enlargement of 
the orbital shell-model space. The level density of states 
0+ in the N = Z nucleus 64 Ge becomes sensitive to the 
inclusion of the next shell (gg/ 2 ) only at the excitation 
energy greater than 14 MeV (see the right graph of Fig. 
0 which means that the region of neutron resonances 


could be reliably evaluated with the more narrow orbital 
space. This case has important ramifications for the as- 
trophysical consideration of the element abundance since 
this nucleus is considered to be a waiting point in the 
r-process of nucleosynthesis. 


V. COHERENT AND INCOHERENT 
INTERACTIONS 

As the whole shell-model Hamiltonian contributes to 
the traces defining the level density, we can explore the 
effects of individual components of the effective interac¬ 
tions including the “incoherent” parts of the full Hamil¬ 
tonian which do not significantly contribute to the for¬ 
mation of the mean field. These parts of the interaction 
determine the finite lifetime of the simple quasiparticle 
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FIG. 12: (Color online.) Interpolation of the single-particle level density parameter a, Eq. (1151) . Left panel: different colors 
present different isotopes or isotones. Right panel: moments-method calculation with interpolation (black circles), fit using 
the experimental data on neutron resonances (Ref. [5ll| . orange diamonds), and fit using experimental low-lying levels (M, 
yellow squares). 


(or collective) modes and their fragmentation in terms of 
genuine complicated eigenstates of exceedingly entangled 
nature. In particular, these collision-like interactions are 
responsible for the formation of chaotic states with high 
information entropy and the process of thermalization 
0. For example, it was shown [~18lIbtll that the exactly 
considered pairing interaction contains some chaotic fea¬ 
tures but they are still not sufficient for establishing the 
complete chaotic picture comparable to the predictions 
of the Gaussian Orthogonal Ensemble. 

Below we show the evolution of the level density as a 
function of the interaction modes included in the full cal¬ 
culation of the moments. Although it is not difficult to 
vary all individual matrix elements, here we use a simpli¬ 
fied approach, presenting the whole two-body interaction 
Hamiltonian in the sd-shell model as consisting of two 
parts with variable intensity, 

H = h + fciH(pairing) + fc 2 H (non — pairing). (9) 

Here h contains the single-particle energies, V (pairing) 
all matrix elements with the pairs of nucleons in the 
channel J n T = 0 + l, while all other shell-model matrix 
elements are attributed to the last term. The numerical 
coefficients k\ and fc 2 are varied giving rise to different 
versions of the shell model; the realistic case emerges at 
k\ = /c 2 = 1. It is easy to understand that in the mo¬ 


ments method (and in the full diagonalization, see 0 ) 
new independent components of the Hamiltonian add in 
quadratures to the final width of the shell-model level 
density. 

The global evolution of the full level density in 28 Si 
as the coefficients ki = /c 2 are varied from 0.1 to their 
realistic values is shown in the left graph of Fig. [5] The 
low residual interaction obviously keeps untouched the 
independent-particle partition structure of the Hilbert 
space that reminds the results of the mean-field combi¬ 
natorics. As the parameters fci,fc 2 grow, the next curves 
show the development of the final picture. The configu¬ 
rational structure is gradually washed out by the resid¬ 
ual interaction leading to the final smooth level density 
discussed earlier. Let us stress that the observed evo¬ 
lution is not a consequence of the superposition of all 
subspaces with different values of J. The individual sub¬ 
space J n = 0 + , the right graph on Fig. 0 demonstrates 
practically the same evolution. 

Fig. [5] describes the situation when the non-pairing 
components of interaction are suppressed, fc 2 = 0.1, 
but the pairing strength evolves. Superposing all val¬ 
ues of J (the left graph of Fig. [9]), we see that the 
smooth Gaussian-like curve is achieved only at the non- 
realistically high pairing strength. Again the picture is 
nearly the same when only the states J = 0 are con- 

















FIG. 13: (Color online.) Level density parameter (1111) fitted 
in the simple Fermi-gas model for sd and pf shell nuclei. The 
empty circles (red color) present the fitting range 1-5 MeV, 
the filled black circles correspond to the fitting range 5-25 
MeV. 



FIG. 14: The relative changes in the level density parameter 
a due to the pairing energy shift <5, Eq. cn shown for sd- 
and p/-shell nuclei. The empty circles present the fit with 
the energy shift 5 = — 1 MeV and the filled circles correspond 
to 5 = 1 MeV. 


sidered, the right graph of Fig. [9] The realistic pairing 
strength, k\ = 1.0, at the absence of non-pairing interac¬ 
tions is still not sufficient for the fully smooth level den¬ 
sity. At low excitation energy < 20 MeV the evolution of 
the level density for J = 0 clearly shows the disappear¬ 
ance of the typical large oscillations with the growth of 
pairing. Here, indeed, the pairing interaction shifts the 
noticeable part of levels to higher energies. If the pairing 
matrix elements are fixed at the empirical value, Fig. 1101 
the large bumps from the original partitions do not ap¬ 
pear but the incoherent interactions very much broaden 
the final result. The generic character of this scenario is 
confirmed by Fig. [II] for 52 Fe. 

VI. THERMODYNAMIC DESCRIPTION AND 
COMPARISON WITH PHENOMENOLOGY 

A. Simple Fermi-gas 

The shell-model Hamiltonian ID starts from non¬ 
interacting particles or quasiparticles, elementary exci¬ 
tations in the mean field of certain symmetry that de¬ 
termines the appropriate quantum numbers of excited 
states. For nuclei, the adequate image is that of the per¬ 
fect two-component Fermi gas. The ground state of the 
system is the filled Fermi sphere, and the excited states 
are described by the particle-hole picture. In the realistic 
many-body physics, this is just an initial step that has to 
be followed by switching on the interaction between par¬ 
ticles and holes. However, already by this mechanism, 
the level density increases exponentially which justifies 
the traditional phenomenological approaches. 

The particle-hole phenomenology uses the steepest de¬ 


scent method to calculate the level density as a function 
of excitation energy E through the Laplace transform of 
the partition function which leads to the standard result 
for one type of particles, 

7 ^^’ {w) 

where the level density parameter is 

7r 2 

a = —v F (11) 

6 

and v F is the density of single-particle states at the Fermi 
surface. The generalization to a proton-neutron system 
leads to a modified expression, 


P(E) 


V^a 2VTe 
12 ( clE ) 5 / 4 


6 1/4 V 2y/aE 

12 (v F Ef/ A 


( 12 ) 


where the parameters a and a now include the total 
single-particle density of states at the Fermi surface, 
i/p = v F (n) + is F (p), and the effective single-particle 
density D = v F /(2 V v f{ji)v f (jp)), correspondingly, see 
Eq. (fill) . The singularities at E —> 0 in densities Util) 
and m show that the statistical method of calculation 
is invalid at too low excitation energy where the num¬ 
ber of nuclear many-body states is small. The notion 
of the level density requires that the excitation energy 
be greater than the average distance between the single¬ 
particle levels, 
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FIG. 15: (Color online.) Thermodynamic temperature as a 
function of excitation energy, Eq. urn calculated in the mo¬ 
ments method. The top four curves (black color) correspond 
to the densities of 28 Si and the bottom four curves (blue color) 
correspond to the densities of 56 Fe. 


FIG. 16: (Color online.) The constant temperature T and 
the “constant” (see Eq. (1191 1 fitted for the energy range 5-15 
MeV. Different sets of nuclei from the pf shell with complete 
and almost complete shells are presented by different colors: 
Z = 20 (red), Z = 21 (green), Z = 22 (blue), and N = 40 
(brown). 


In thermodynamic language, the nuclear temperature t 
is introduced through the Fermi-gas formula for the ex¬ 
citation energy, 

E = at 2 . (14) 

In general, for the low excitation energy region it is al¬ 
ways recommended to use directly the available experi¬ 
mental information. 

If one tries to compare the thermodynamic level den¬ 
sity of the Fermi-gas type with experimental data, it is 
hard to expect the numerical agreement of the level den¬ 
sity parameter m with that required by data even if 
the exponential growth of the total level density takes 
place. As mentioned in the Introduction, when the level 
density grows, the residual interactions lead to multiple 
avoided crossings and mixing of many-body levels with 
the same exact quantum numbers. This process of chao- 
tization evolves the level network considered as a function 
of the interaction strength close to the aperiodic crystal 
with a small average spacing. The whole set of stationary 
states is becoming locally close to the predictions of the 
Gaussian Orthogonal Ensemble of random matrices. The 
wave functions here are quite complicated superpositions 
of very many particle-hole states. The energy behavior of 
the level density is now close to a Gaussian Q with the 
total width that is given by adding in quadratures the 
initial width due to the mean-field quasiparticles and the 
dispersion of the off-diagonal matrix elements of residual 
interactions. 

Fig. IT?1 shows the parameters a (the top panels) and 
constant (the bottom panels) found from the shell-model 


calculation of the total level density fitted by 

In [p(E,M = 0)] = 2 \faE — ^ In E + const. (15) 

Both the parameter a and the constant show the change 
clearly correlated with the microscopic filling of the nu¬ 
clear shells. The parameter a reveals the maximum in 
the middle of the shell occupation as it should have been 
expected from the construction of the model. The em¬ 
pirical estimates for the same nuclei are available from 
the level density at the neutron resonances energy |5lj . 
which do not show considerable shell effects, and by ex¬ 
trapolation from low-lying levels (52j where one can see 
very weak shell effects in the region of the mass num¬ 
ber around A ft: 50. The constant in Eq. 113 is small 
but also shows in some cases the shell-model dependence 
with a minimum in counter-phase with the parameter a. 
Fig. M shows the dependence of the level density pa¬ 
rameter a on the energy range for which the fitting was 
performed. The empty circles (red color) present the low 
energy fit, around 1-5 MeV, while the filled circles (black 
color) present the standard 5-25 MeV energy range fit. 
We can see that at low energies the fitted level density pa¬ 
rameter a is slightly larger, but still it is not large enough 
to be compared with the empirical estimates j5l], [53] ■ 

B. Back-shifted Fermi-gas model 

In this paper we do not discuss in many details the 
nuclear pairing correlations and its importance for the 
nuclear level densities leaving this for the future consid¬ 
eration. One of the standard phenomenological way to 




















FIG. 17: Comparison between global level densities for 28 Si 
calculated with the moments method (the solid curve) and 
using the constant temperature for 5-15 MeV (the dotted 
curve) and for 5-25 MeV (the dashed curve) energy range; 
fits according to Eq. m- The inset presents the enhanced 
low-energy region. 



FIG. 18: (Color online.) Constant temperature fitted for sd 
and pf shell nuclei. The diamonds (red color) present the 
fitting range 1-5 MeV, the circles (black color) correspond to 
5-15 MeV and the squares (blue color) correspond to 5-25 
MeV fitting ranges. 


account for pairing correlation is to use the back-shifted 
Fermi-gas formula (BSFG) fl8l. l53i| . 

p t E)= _ El _ 2 (16) 

Pbsfg(u) i2 [a(E — < 5)] 5 / 4 6 ’ 1 j 

where the excitation energy is shifted by the pairing en¬ 
ergy parameter < 5 . This introduces a new free parameter 
that should be fitted alongside with the level density pa¬ 
rameter a. The pairing energy parameter 8 is in general 
different for even-even, odd-odd, and even-odd nuclei due 
to the formation of Cooper pairs, and does not necessarily 
coincide with the corresponding pairing gap parameters 
A or 2A. 

Fig. [14] shows the sensitivity of the level density pa¬ 
rameter a to inclusion of the energy shift 8 . The quantity 
plotted along y- axis is the relative change in a if we shift 
the excitation energy by 8 , 

100% x a ~ aBSFG , (17) 

a 

where clbsfg is the level density parameter fitted using 
Ea. dTfil) with a fixed value of the pairing energy 8. We 
can see that the a parameter varies only in a ±6% range 
when the shift changes from 8 = — 1 MeV (presented by 
the empty circles) to <5 = 1 MeV (filled circles). 


C. Constant temperature model 

The model of the energy dependence of the level den¬ 
sity different from the Fer mi-g as phenomenology (1121) be¬ 
ing suggested long ago [18|,|5S| gradually becomes popular 
among practitioners. It is assumed that the level density, 
at least up to 10 MeV excitation energy, and maybe even 


higher Ihil lbC] , can be described by the constant tempera¬ 
ture T. This temperature is the single parameter defined, 
in the simplest version, according to the thermodynamics 
as 


din p(E) 
dE 


(18) 


The philosophy behind this approach is usually explained 
[57j in terms of the first-order phase transition that goes 
through the latent heat at fixed temperature. Although 
typically this assumes the melting of the Cooper pairs but 
in fact one can also talk about other types of correlated 
structures which are undergoing something similar to the 
liquid-gas phase transition or even the first stage on the 
road to multifragmentation. In a more detailed descrip¬ 
tion, the effective temperature parameter can be different 
for the classes of states with different quantum numbers, 
although such a generalization does not look well from 
the viewpoint of the thermal equilibrium between vari¬ 
ous degrees of freedom. Such an effective temperature 
parameter (plus a corresponding constant) could be fit¬ 
ted in a reasonable energy range to represent the partial 
(with certain spin and parity) or total nuclear level den¬ 
sities as 


E 

In [p(E, J)] = — + constant, 
Tj 

E 

ln[p(.E)] = — + constant. 


(19) 


It is immediately clear that globally the constant tem¬ 
perature model cannot be compatible with our shell- 
model calculations. In the truncated orbital space the 
global level density will always look as a Gaussian with 





























FIG. 19: (Color online.) Logarithm of level density, In [p(E,M)\, versus M 2 , Eq. (1231) . calculated for 28 Si, 52 Fe, 44 Ca, and 
64 Cr. Different colors present different excitation energies: 5 MeV (black), 10 MeV (red), 15 MeV (green), 20 MeV (blue), 
and 25 MeV (brown). Dots correspond to the moments method calculations, the solid lines are linear interpolations, see the text. 


the effective temperature (fT51) 


changing with energy from positive to negative values 
on different sides of the centroid energy E c . Here erg 
is the effective width of the Gaussian that reflects the 
summed contribution of all components of interactions. 
However, at relatively low energy the level density can 
grow approximately exponentially effectively resulting in 
an approximately constant temperature. 

Here we give a couple of examples showing the ex¬ 
ponential fit to the level density of different j-classes 
and global level density. Fig. [15] shows how the actual 
temperature (see Eq. (fI51) l calculated by the moments 
method for 28 Si and 56 Fe depends on the excitation en¬ 
ergy. The temperatures calculated for certain spins J and 
for all spins (the total density) are not constants, they 
increase with the excitation energy as suggested in Eq. 

The corresponding constant temperature fit of Eq. 
(1191) performed for the sd and pf nuclei is presented in 


Fig. [TC] The effective constant temperatures in the figure 
depend on the range of the excitation energies where the 
fit was performed, — the greater the excitation energy 
the higher the effective constant temperature. 

We note that the majority of the fitted temperatures in 
Figs. flSlandflblare concentrated in small regions near 2-5 
MeV, while there are some exceptional cases of unreason¬ 
ably high temperatures of 10-20 MeV that correspond to 
the nuclei with complete or almost complete shells for one 
or two sorts of nucleons when the Fermi gas approxima¬ 
tion is obviously invalid. Finally, Fig. [T71 shows how good 
the constant temperature approximation is. The dotted 
curve (fitting energy range is 5-15 MeV) and the dashed 
curve (the range is 5-25 MeV) present the correspond¬ 
ing constant-temperature level densities for 28 Si. We can 
see that these densities work pretty well inside the en¬ 
ergy interval where they were fitted compared to the “ex¬ 
act” density calculated within the moments method. As 
the excitation energy increases the constant-temperature 
densities stop working going too high very fast. Fig. 
m shows the fitted temperatures when different energy 
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FIG. 20: (Color online.) Interpolation of the spin cut-off parameter cr 2 , Eq. (1241) . Parameters a. and (5 are shown for different 
sd and pf nuclei. Left panel: different colors present different isotopes or isotones. Right panel: moments-method calculation 
with interpolation (black circles), statistical calculation (orange diamonds), rigid sphere approximation (yellow squares). 


ranges are used: red diamonds present 1-5 MeV fitting 
range, black circles and blue squares present 5-15 MeV 
and 5-25 MeV energy ranges correspondingly. It is nat¬ 
ural that at higher excitation energies the effective tem¬ 
perature increases thus reducing the rate at which the 
density grows, see Eqs. and 420]). 

VII. SPIN CUT-OFF PARAMETER 

The distribution of levels with certain values of global 
nuclear constants of motion (angular momentum, parity, 
isospin) is of special interests in all applications. Our 
method directly supplies the required information for ev¬ 
ery set of those exact quantum numbers. The standard 
phenomenological approach to the angular momentum 
dependence of the level density assumes the random an¬ 
gular momentum coupling as a diffusion process in the 
space of projections M. The total projection results from 
the random walk, and the fraction of states with a given 
projection M is Gaussian, 

P(E,M ) _ 1 —M 2 /2a 2 / 91 n 

p(E) 1 J 

Then, as it was mentioned earlier, the density of states 
with a given value of J is just a difference 

pj(E) = p(E,M = J)-p(E,M = J+ 1). (22) 


Assuming this random angular momentum coupling we 
expect the linear M 2 -dependence of the logarithm of the 
level density, 

M 2 

In [p(E, M)} = In \p(E, M = 0)] - ^. (23) 

The top two graphs of Fig. H9lshow the M 2 -dependences 
of In [p(E, M )] at different energies for 28 Si and 52 Fe. The 
lines correspond to the best linear interpolation of this 
logarithm. In these examples we see a very good linear 
behavior and the smooth dependence on excitation en¬ 
ergy. This can be interpreted as an evidence for random 
coupling of angular momenta of individual particles. 

The situation is different in two next examples, the 
bottom graphs in Fig. [T9l 44 Ca and 64 Cr. Here we do 
not have a regular energy dependence and, therefore a 
clearly defined parameter cr. This is what could be ex¬ 
pected from physical arguments. The whole idea of the 
Gaussian random walk in the angular momentum space 
breaks down here because of the isospin limitations. The 
nucleus 44 Ca in the shell-model description has only four 
identical / 7 / 2 neutrons which allow for the isospin T = 2 
only and for the interaction in the particle-particle chan¬ 
nel with the total isospin T = 1. Therefore many val¬ 
ues of the total spin are forbidden. The second nucleus, 
6 4 Cr 40 , in the p/-shell model has only four valence pro¬ 
tons with the same limitations of the angular momentum 
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coupling. In both cases, it is hard to expect the require¬ 
ments of the random spin coupling to be satisfied. 

The value of the spin cut-off parameter a can be ex¬ 
tracted from the curves as in Fig. ITTjl where we observe 
a good linear behavior. Using the thermodynamic lan¬ 
guage, we expect this parameter to be proportional to 
temperature, or, for Fermi-gas, to the square root of en¬ 
ergy. The corresponding parameterization can be taken 
as 

a 2 = ol\/~E (1 + (3E). (24) 

The coefficient a can be taken 0 from Fermi gas sta¬ 
tistical mechanics as oc ppT(M 2 ), where vp is the single¬ 
particle level density at the Fermi surface, or assuming 
the angular momentum corresponding in average to the 
rigid-body rotation with the moment of inertia oc 7L4 5 / 3 . 
The results of the shell-model calculations are shown in 
Fig. [20l We define the parameters a and /3, Eq. (l24l) . 
from the energy region 5-25 MeV. Comparing two groups 
of nuclei, sd-shell and pf-shell, we indeed see the average 
growth of the spin-cut off parameter for two represen¬ 
tative groups as proportional to A 5 / 3 . It is impossible 
here to make a selection between the statistical estimate 
of the spin cut-off parameter and the estimate from the 
moment of inertia as both of them, being too crude to 
reflect shell efects inside each group which are certainly 
present, do not agree with the A 5 / 3 estimate and require 
a more detailed analysis. The constant /3 from Eq. (l24l) 
is small but, at least for the pf- nuclei, may also contain 
some shell effects. 


VIII. CONCLUSION 

In this article we collected, explained and over viewed 
the first results of the improved method for statistical 
calculation of the nuclear level density for a given shell- 
model Hamiltonian. The method is physically based on 
the chaotization of the intrinsic dynamics by the inter¬ 
particle interactions. In practice, one needs to calculate 
only the lowest moments of the Hamiltonian partitioned 
in terms of mean-field configurations. The first two mo¬ 
ments turn out to be sufficient for the full agreement of 
the found level density with the result of the exact di- 
agonalization as it was checked by the cases when such 
full diagonalization was technically possible. The serious 
improvements compared to the previous attempts in the 
same direction include the use of the finite-range Gaus¬ 
sian distributions and of the recurrence relation for elim¬ 
inating the spurious states. We did not discuss the de¬ 
termination of the ground state energy that is necessary 
for the appropriate positioning of the level density. There 
are special methods for doing this, including the exponen¬ 
tial extrapolation also based on the chaotic properties of 
remote highly excited states [13 . The shell-model level 
density can be calculated in any specific class of global 
constants of motion (proton and neutron numbers, to¬ 
tal spin, parity and isospin) as a function of excitation 


energy. This is essentially what is needed for practical 
applications to nuclear reactions including those in as¬ 
trophysics. 

The main conclusion that can be drawn from this ex¬ 
perience is that the shell-model level density that re¬ 
sults from the statistical calculation is a smooth function 
of excitation energy in all classes of quantum numbers 
containing a considerable number of states allowed by 
the truncation of the orbital space (of course, the state 
with the maximum possible total spin is frequently just 
unique). Giant oscillations of the level density predicted 
by the calculations based on the mean-held combinatorics 
are almost completely erased by the presence of incoher¬ 
ent collision-like interactions which usually remain out¬ 
side of the mean-held models or parametrizations with 
the so-called collective enhancement. Taking into ac¬ 
count all components of residual interactions, coherent 
(such as pairing) and incoherent, is necessary for the ad¬ 
equate description. 

The comparison with phenomenological Fermi-gas ap¬ 
proaches, including the models with constant tempera¬ 
ture, shows that, being less theoretically justihed than 
the full direct calculation, in many cases they are nev¬ 
ertheless quite reasonable for practical use. The calcu¬ 
lation of the spin dependence of the level density and of 
the relevant spin cut-off parameter is more sensitive to 
assumptions, and there are cases when it is not in good 
agreement with exact results. The model of constant 
temperature, in our opinion, can be applied at relatively 
low excitation energy but, most probably, it reflects the 
general process of chaotization the dynamics rather than 
just breaking of Cooper pairs. Certainly, the accumu¬ 
lation of experimental data and new applications of the 
statistical method are necessary for better understanding 
the underlying physics. 

The whole approach unavoidably suffers from the gen¬ 
eral problems of the shell model. It is possible to be¬ 
lieve that the results will not be sensitive to the specific 
versions of the shell-model Hamiltonian as soon as this 
choice agrees well with the low-lying spectroscopy. How¬ 
ever, the space truncation provides a natural limitation 
for the applications of all such methods. The space can 
be expanded (and many-body residual interactions can 
be included) paying the price of longer computational 
time that can be cut off by the parallelization. But any¬ 
way for any choice of finite space there is a natural limit 
of applicability. Luckily enough, it seems that this limita¬ 
tion is not essential for many astrophysical applications. 
More theoretical work is necessary for understanding the 
resonance density for the states deeply in the continuum 
which again might not be critical for a typical stellar tem¬ 
perature when the resonance states under consideration 
still are quite narrow. 

The whole development of the method was done in 
collaboration with M. Horoi. The discussions with B.A. 
Brown are acknowledged. The work on level density was 
supported by the NSF grants PHY-1068217 and PHY- 
1404442. 
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